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We prove that it is possible to freeze a light pulse (i.e., to bring it to a full stop) or even to 
make its group velocity negative in a coherently driven Doppler broadened atomic medium via 
electromagnetically induced transparency (EIT). This remarkable phenomenon of the ultra-slow 
EIT polariton is based on the spatial dispersion of the refraction index n{ui, k), i.e., its wavenumber 
dependence, which is due to atomic motion and provides a negative contribution to the group 
velocity. This is related to, but qualitatively different from, the recently observed light slowing 
caused by large temporal (frequency) dispersion. 



PACS numbers 42.50.Gy, 42.55.-f, 42.65.Wi 



Slow group velocity in coherently driven media [Q has 
been shown to provide new regimes of nonlinear interac- 
tion with highly increased efficiency even for very weak 
light fields, high precision spectroscopy and magnetome- 
try H . It has been demonstrated [Q^Ql that EIT is ac- 
companied by large frequency dispersion, \ujdn/dLu\ 1, 
and can slow the group velocity down to 10 - 10 2 m/s. 

In this paper we show that, using spatial dispersion 
due to atomic motion, it is possible to freeze the light 



Using the convolution theorem to write the RHS as 
J dkduxi^! k)£{ui, k) exp i(Qt — kz), expanding the sus- 
ceptibility to the first order in k, Co, noting that k and Co 
under the integral may be written in terms of d/dz and 
d/dt acting on £(t, z) and rearranging terms we have 



0, or even to make its group velocity opposite to the which im P lies the field equation with v g given by Eq. fllj) , 



wavevector, v g < (see Eq. ([[])). We consider two differ- 
ent types of atomic media: (i) atomic beam or uniformly 
moving sample, and (ii) hot gas in a stationary cell. 

Freezing of light in a stationary cell via a hot gas is es- 
pecially intriguing (Fig. 1). The idea is to tune the driv- 
ing field to resonance with the velocity group of atoms 
that moves in the direction opposite to the light pulse 
with velocity equal to the light group velocity that would 
be supported by this group of atoms if they were at rest. 

The main result of the present paper is contained in 
Fig. 2 which shows that v g can be zero, for a pulse in a 
hot gas, when the drive detuning Alu c i is properly chosen. 

As is well-known, in a medium possessing both tem- 
poral and spatial dispersion of the refraction index, 
n(uj, k) = 1 + 2irx(u, k), the group velocity of light con- 
tains two contributions, 
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Eq. (Q) is an immediate result of differentiating the 
dispersion equation kc = u>n(to,k), i.e., c = v g (n + 
ujdn/duj) + uodn/dk. The meaning of Eq. (|J) becomes 
clear if one turns to the equation for a field amplitude 

(^ c ~g — I" ~g^j ^ = 2tt{uj J dz'dt'x(t — t' , z — z')£(t' , z'). 



The first term in Eq. (0), v g = Re[c/(n + ujdn/ du)}, 
is due to frequency dispersion, and was discussed 
in recent papers (lJ-H]. The second term, v s = 
Re{(udn/dk)/(n + udn/duj)], is due to the effect of spa- 
tial dispersion, i.e., nonlocal response of the medium to 
a probe field. We study dilute systems where the suscep- 
tibility is small, |x(w, fc)| <§; 1, but v g <C c, as it is for 
all the EIT experiments carried out so far. As usual, we 
consider real- valued group velocities under the condition 
that imaginary part of duj/dk is negligible. Otherwise 
group velocity looses its simple kinematic meaning and 
strong absorption governs or prevents propagation of the 
light pulse through the medium. The latter is the rea- 
son why the resonant interaction of light with a two-level 
medium never results in an ultra-slow polariton. 

A mono-velocity atomic beam or uniformly mov- 
ing sample corresponds to the simple case of spatial 
dispersion, so-called drift dispersion. In the co-moving 
frame atoms are at rest, there is no spatial dispersion 
and the group velocity is given by the first term of Eq. (1) 
alone, v g . The Galilean transformation to the laboratory 
frame, k = k, lu = uj — kv, where v is the atomic velocity, 
yields the group velocity v g — Re(duj / dk) = v g — v. 
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Eq. (1) yields the same result, since the susceptibility 
i|l depends only on the combination ui + kv, 

ip 2 ab N n ab T cb 



X«(w, k) = x(w + feu) 



- ^ 2 n CQ /r* c 



r a 6r C f, + Q 2 



(2) 



Here n ab = p aa - p bb , n ca = p cc - p aa , p u is the pop- 
ulation of the ith level, 7 and 7 C (, are the relaxation 
rates of excited state and c — b coherence respectively 
(7 7cb); u ab and u> cb are the frequencies of the op- 
tical and low frequency transitions (uj ab 3> w cb ); w<j, 
fed and u>, fe are the frequency and wavenumber of the 
driving and probe fields respectively, N is the atomic 
density, Q — \ii ac E^\/2% is the Rabi frequency of drive 
field (1/2) Ed ex.p(iu)dt — ik^z) + c.c, p ac and p a b are the 
dipole moments of a — c and a — b transitions respec- 
tively, T ac = 7 + i(Atv d + k d v), T ab — 7 + i(Aw + kv), 
r c b = 7cb + »(Aw - Aw d -I- Afeu), Auj d = u d - Lu ac , 
Auj = uj — u> ab , kd — uj d /c, k = k d + u cb /c + Ak. We use 
a standard model with incoherent pump and loss rates 
(v c = r b — 7c&/2), assuming time of flight broadening 
of b — c transition (Fig. 1), so that in the absence of 
fields p cc = p bb = 1/2. According to Eqs. (1), (2), we 



again obtain v g 



v. The physical reason for this 



drifting of the pulse is that the field is basically "seized" 
by the atoms in the form of atomic polarization. 

An important question is how to input the light pulse 
into the gas. There are different possibilities. One exam- 
ple uses a grid mirror that has the grid stripes of small 
area, so that atoms can freely fly through the mirror, 
and small spacing between the grid stripes as compared 
to the wavelength of light to provide efficient reflection, 
as in Fig. lb. If atoms are at rest, the light would prop- 
agate in the forward direction. However, if the velocity 
of atoms is equal to (or larger than) v g , one should see a 
frozen (or backward) pulse. 

Depending on the mechanism of pulse input into the 
medium, one should look for the solution of the problem 
with initial (time), boundary (space), or mixed (time- 
space) conditions. In the case of the initial value prob- 
lem, we solve the dispersion equation for u> = to(k), 
Fig. 3a. Galilean transformation ensures the same EIT 
half width, AksiT = Q 2 /jv g , as for the atoms at rest 
since Im[u>(k)] = Im[u>(k)]. In the case of the bound- 
ary value problem, we find k = k(u>). The result shows 
narrowing of the EIT dip proportional to the kinematic 
factor a = (v g — v)/v g . Indeed, in the accompanying 
frame the dispersion relation near EIT resonance can 
be decomposed in a form of a quadratic polynomial, 
Afe = Ak Q - i[n + (,(ACj - Aui d ) 2 } + (Aw - Auj d )/v g . Its 
Galilean transformation to the laboratory frame yields 



Afe = Afe + - 
a 
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(3) 



where Slo — Auj — Au>d- Coefficients in Eq. (3) can be 
easily deduced using Eq. (2). For example, for the case 



of one-photon resonance Aui d = at fl 2 3> 7ch7, we have 
v g = hQ 2 /2Trp ab kdN, Ak — 0, residual absorption co- 
efficient at the center of EIT dip is kq = ^cb/vg, and a 
coefficient determining the parabolic profile of absorption 
in the EIT dip is £ = j/fl 2 v g . This approximation is valid 
if residual absorption is small, k £ <C (1 — v/v g ) 2 /v 2 . Ab- 
sorption increases twice as much as EIT minimum value 
at detuning Slue it = \Vg — v \^yldJl / v g that is much 
less than the EIT half width Alueit — AkEiT\v g — v\. 

Eq. (3) shows that the absorption coefficient Im k is 
increased and sharpened by a factor (v g — v) jv g as com- 
pared to that in the co-moving frame. Since the spec- 
trum of the pulse cannot be transformed on the station- 
ary boundary, only those spectral components that are 
within the sharpened EIT dip penetrate deep into the 
medium. For drift velocity v > v g , the backward EIT 
polariton can be excited from inside a cell (Fig. lb). 

In the case of an atomic beam with a moving boundary 
(or moving sample), i.e., for the mixed boundary-initial 
value problem, the spectrum (inverse duration) of the 
pulse shrinks at the moving boundary exactly in the same 
way as the EIT width in Eq. (3), Aw = Auj(v g — v)/v g . 
This is not a coincidence, but is necessary for consistency 
of viewing of the same process from different frames. The 
pulse within the EIT dip decays in time with the same 
rate independently of whether it propagates through 
atoms at rest or through a beam, since this decay is pre- 
determined by atomic relaxation 7 C &,7- 

Atoms with a thermal velocity distribution. Let 
us consider a stationary cell of hot atoms. If the intensity 
of the drive is strong enough to provide EIT for the res- 
onant group of atoms (see Fig. 4) but at the same time 
weak enough to avoid an interaction with off-resonant 
atoms, moving with "wrong" velocities, it is mainly this 
drifting beam that would support the ultra-slow EIT po- 
lariton with zero or even negative group velocity. 

To prove this we calculate the dispersion law u(k) for 
the EIT polariton in a hot gas in a cell at rest. The sus- 
ceptibility is given by an average of the beam susceptibil- 
ity over a velocity distribution F(v) of atoms in a gas with 
thermal velocity vt, x(w,fc) = dvF(v)xv{u, k). In- 
stead of the Maxwellian thermal distribution we can use 
Lorentzian, F(v) = vt/[t^(v t + v 2 )], since the far-off- 
resonant tails are not important. This allows us to obtain 
simple analytical results because an integration over ve- 
locities is reduced to a sum of a few residues in the simple 
poles, v = Vj. Only those poles count that lay in the lower 
half complex u-planc in the formal limit of infinitely large 
growth rate Imu — > — 00. For a positive wavenumber de- 
tuning, Afe > 0, there are two such poles. One originates 
from Lorentzian, V\ = —ivx, and the other from the ve- 
locity dependent populations, V2 = —(ijG + Aujd)/kd- 

1/2 

Here 7G = j (l + Q 2 / (^cbl)) determines the veloc- 
ity width of an effective drifting beam of atoms that 
are driven by an external field into a coherent "dark" 
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state an d> hence, responsible for the ultra-slow 

EIT polariton (see Fig. 4). For Ak < 0, there is an 
additional pole, v% oc 1/Ak, originated from resonance 
FabFcb + Q 2 — in Eq. (2). However, near EIT reso- 
nance, i.e., for small detuning Ak, it enters the lower half 
plane from infinity, v% — ► — too, so that its contribution 
is negligible if N < fc| (job/l) \/ k d v T /tt. 

Calculation of the residues at poles V\ and V2 yields 



Auj = Au d - v d Ak + - 
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shown in Fig. 3b. The EIT half width is Ak' EIT = jk/v' g - 
For small detuning |Afc| <C Ak' EITl Eq. (g) yields lin- 
ear dispersion and parabolic decay profile, Auj — Auj d + 

~ 2 

Ak{v' g — v d ) + i"f c b + iAk 2 v' /"/k- Decay increases twice 
as much as EIT minimum value, ImAuj — 2"f c t, at very 



, (4) 



small detuning Sk^ 



EIT 



y^fcVu/v'g < Ak' L 



EIT- 



where m = [RiT£) - r$(l + 2 7 i? 1 / 7cb )]/[l + 7 2 (G 2 



The 

group velocity describes pulse kinematics if duj/dk has 
negligible imaginary part, i.e., near the center of the EIT 



dip where |Afc| < 



Vd\lk/v' q . The last inequality 



l)i?i/f2 2 ], 772 = k d v T R 2 [fl 2 /(G - 1) - r^7]/7 c b7G, does not mean that the pulse cannot be stopped. It just 
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7 + kvT + iAuj , r 
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7(1 + Gk/k d ) + i(Au 



kAui d /k d ), ric = 7 + k d v T + iAuj d , T 



(i) 



leb 



\Ak\v T + i(Auj-Auj d ) ,T<§ = lcb + \Ak\jG/k d +i(Auj- 
Au d - AkAuj d /k d ), i?! = n 2 /[~f 2 + (Au d -ik d v T ) 2 }, 
Ri = Q 2 /[{k d v T ) 2 + (Aiu d + i 7 G) 2 ]. 

The susceptibility (4) of a hot gas looks like the suscep- 
tibility of a medium consisting of just two mono- velocity 
components: (i) broad background with velocity v = 
and linewidth 7 + k d vx, and (ii) a drifting beam with 
velocity v d — —Au d /k d and power broadened linewidth 
7(1 + G) (see Fig. 4). This interpretation becomes very 
accurate near EIT dip, \Aui — Aui d \ <§; 7G, at the con- 
ditions necessary for the existence of freezing ultra-slow 
EIT polariton: a) low-frequency coherence decay is much 
slower than optical decay (7^ <§; 7); b) drifting beam 
width is less than Doppler broadening ( 7 G <§; k d vr)\ c) 
detuning of driving and probe fields from one-photon res- 
onance is large enough (lAw^l ^> 7G) while two-photon 
resonance is maintained. Then, for the ultra-slow EIT 
polariton, the susceptibility is approximated as 



(Ml 

hjG 



n 2 



j{l + G){UJ - LU k ) 



(5) 



if we keep only resonant w-dependence in denomina- 
tors setting everywhere else Aw = Au d . Here N' = 
N"fGk d VT/[(k d vr) 2 + Auj 2 ] <C N is the density of atoms 
in the drifting beam. The resonant denominator, where 
uj k = uj ab + Auj d k/k d + ij k , 7^ = 7 c6 + fi 2 /7(l + G) + 
\Ak\jGfkd, comes from the factor Q 2 + r^r^ in 
Eq. (4). Thus, we explicitly find the frequency and the 
decay (u>k, 7&) of the EIT exciton coupling of which to 
the probe field produces the ultra-slow polariton. 

For the boundary value problem, Eq. (5) yields a dis- 
persion that is similar to that for the mono-velocity 
beam (3) with parameters v = v d , v' g = [(k d VT) 2 + 
Au#fi 2 ty[/&iV7(l + G)k 2 v T ], «o = jcb/%, $ = lh k v' g . 

For the initial value problem, from the dispersion equa- 
tion fcc = oj(l + 2nx) and Eq. (5), we find dispersion law 



means that when the pulse is frozen, v„ 



v d 



0, 



its evolution is governed by the dispersion of absorption. 

Fig. 3 clearly shows that the ultra-slow EIT polariton 
in a hot gas is similar to that in a mono-velocity beam, 
since detuning of driving field picks a beam with velocity 
v d = —Auj d /k d . However, effective density of atoms sup- 
porting EIT polariton N' and EIT width Ak' EIT = "/ k /v' g 
in a hot gas are different because of factors 7G and F(v). 
As a result, the group velocity at the EIT resonance, ac- 
cording to Eq. (6), in terms of a critical density is 



0N C , 



NF(v d ) 



Vd, N c 



(7) 



where (3 — max[v d F(v d )]. For Lorentzian F(v d ), we have 
13 = l/2vr, and v g = (v d - v d ly )(v d - v { d 2) )N cr /2Nv T 

is a quadratic polynomial over v d , i.e., the group ve- 

(1 2) 

locity is zero for drive detunings v d — vt[N/N ct ± 
■J (N/N cr ) 2 — 1] and negative between them for den- 
sity higher than the critical value, N > N cr , as is 
shown in Fig. 2. To achieve minimal group velocity, 
mini; 5 = -{v T N/2N cr )[l - (iV cr /iV) 2 ], one has to tune 
at v d = VTN/N cr . The condition to freeze or reverse the 
light (v g < 0) means that the group velocity supported by 
the drifting beam with the density N' — ■nNF(v d )"fG jk d 
should be equal to or less than the velocity of atoms in 
the beam, i.e., v g = v g N'/N < v d . If we compare a 
mono-velocity beam with a hot gas at v d = v and the 
same N' as the total density N in a beam to provide the 
same group velocity, v g = v' g , we find that the EIT width 
and the residual decay in a hot gas are G ~ ^/V^wy 
times less than in a beam. To minimize N cr the drive 
intensity should be as low as possible to decrease v' g due 
to power broadening effect and to avoid EIT contribution 
from the atoms with "wrong" (positive) velocities. That 
is, the drive intensity should be just above a threshold 
of the EIT effect at resonance, il 2 > J c bl- Under real- 
istic parameters relevant to the experiments with 87 Rb 
vapor Q and chozen in Figs. 2-5, the critical density is 
N cr - 10 11 cm" 3 . 

Absorption or time variation of the drive field results in 
a spatial or time dependence of the group velocity in the 
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cell. This allows us to control input and parameters of 
the pulse in the cell. According to geometrical optics, the 
parameters of the EIT polariton adiabatically follow the 
local properties of the driven atoms. Fig. 5 demonstrates 
how the ultra-slow pulse decelerates up to the point v g — 
where it becomes frozen. 

The important conclusion is that the drifting beam 
provides large enough drift spatial dispersion dn/dk (see 
Eq.(l)) to ensure v g < 0. Although the density of drifting 
atoms is small N' <C N, their resonant contribution dom- 
inates. This allows us to make the group velocity zero or 
even negative H. To observe freezing or backward light 
one can look, e.g., for a scattering, luminescence, delay, 
or enhanced nonlinear mixing caused by ultra-slow pulse. 
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FIG. 2. Ultra-slow and negative group velocity of EIT 
polariton vs detuning of drive laser; fl = 0.257, kdVr = IOO7, 
lcb = O.OOI7, (a) N = 0.6N cr ; (b) N = N cr ; (c) N = 1.5iV cr . 



FIG. 3. Dispersion, ReAuj = Re[iv — ui a b], an< i decay, 
ImAuj, spectra of the ultra-slow EIT polariton according to 
numerical solution of the dispersion equation for: (a) atomic 
beam (N = l.lN cr nF(v d )-yG/k d ) with susceptibility (2); (b) 
stationary cell of hot gas (N = l.ljV cr ) with exact susceptibil- 
ity (4); Q = 0.257, v = v d = vt, k d v T = IOO7, -y cb = O.OOI7. 



FIG. 4. The velocity distribution of atoms in a cell (solid 
line) . Effective drifting beam (dotted) selected by drive laser. 



FIG. 5. Kinematics of the deceleration of the ultra-slow 
pulse to the point of freezing (v g = 0) along a cell with de- 
creasing group velocity v g (z). Positions of pulse are shown 
at subsequent moments of time t = mr (r = 3L/2v g (0), 
m=0,l,2,3). v g (z) is calculated numerically according to de- 
creasing drive intensity found from the wave equation for the 
same parameters as in Fig. 3(b), L = 10 cm. 



FIG. 1. (a) Three-level atomic A-system. (b) Geometry 
of ultra-slow EIT pulse propagation in the gas of atoms. 
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